Fast optimal transition between two equilibrium states 
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We demonstrate a technique based on invariants of motion for a time-dependent Hamiltonian, 
allowing a fast transition to a final state identical in theory to that obtained through a perfectly 
adiabatic transformation. This method is experimentally applied to the fast decompression of an 
ultracold cloud of Rubidium 87 atoms held in a harmonic magnetic trap, in the presence of gravity. 
We are able to decompress the trap by a factor of 15 within 35 ms with a strong suppression of the 
sloshing and breathing modes induced by the large vertical displacement and curvature reduction 
of the trap. When compared to a standard linear decompression, we achieve a gain of a factor of 37 
on the transition time. 
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I. INTRODUCTION 

The controlled manipulation of quantum states is cen- 
tral to many areas of physics such as quantum informa- 
tion processing [l|, , design of pulses for nuclear mag- 
netic resonance imaging , atomic gas cooling |6'| and 
transport [7j, or ion manipulation [8]. The paradigm of 
adiabatic transformations, in which the Hamiltonian pa- 
rameters are changed infinitesimally slowly with time [9| , 
is often used to drive a system from a given quantum 
state to another. However, the urge to shorten the du- 
ration of the experiments has driven the search for fast 
optimal non-adiabatic strategies [Tol - [l3| , with a mini- 
mal amount of extra energy supplied to the system. 

In the particular field of cold atoms, time-dependent 
potentials are becoming increasingly used. Examples in- 
clude the transport of cold atomic samples over various 
distances 0, [U [ll] , or the production of very low tem- 
peratures using trap decompression Q . To minimize the 
energy imparted to the atoms, most of these experiments 
were performed in the adiabatic regime where the process 
duration was much longer than the oscillation period in 
the potential, yielding times in the few seconds range or 
longer. Achieving a faster transfer with a limited heat- 
ing motivated experimentalists to employ various non- 
adiabatic procedures [ll|, • Recently, a method based 
on invariants of motion was proposed for the decompres- 
sion of harmonic traps which was argued to g ive access 
to shorter times than "bang-bang" control [12| , provided 
that negative curvatures could be transiently applied . 

We present in this article an experimental demonstra- 
tion of shortcuts to adiabaticity based on this method, 
which we employ to decompress a cloud of magnetically 
trapped ®''Rb atoms. Because of gravity, the position 
of the trap center shifts vertically, which induces slosh- 
ing modes of the trapped atoms. At the same time, a 
breathing mode is excited by the reduction of the trap 
frequency. We thus generalize the approach of Ref. [11] to 



the case of a time-dependent harmonic plus constant lin- 
ear potential to account for gravity. We derive a trap fre- 
quency trajectory ujz{t) which yields a final state identi- 
cal to that obtained through a purely adiabatic transfor- 
mation (hence the "optimal" transition), but in a much 
shorter time. We experimentally implement this trajec- 
tory to perform a vertical trap decompression by a fac- 
tor of 15 within 35 ms (corresponding to roughly half 
the decompressed trap oscillation period), with a strong 
suppression of the cloud's center of mass and size oscil- 
lations. 



II. OPTIMAL TRAJECTORY 
DETERMINATION 

We start our theoretical approach by considering a 
time-dependent harmonic oscillator in the presence of 
gravity 
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with initial and final angular frequencies uJz{0) = wqz and 
ujz{tf) — ujfz, respectively. The objective is to engineer 
a trajectory U!z{t) between these two values so that if we 
start with an initial state at equilibrium at temperature 
To, this state is mapped to a final equilibrium state at 
temperature Tf = with 7^ = wq^/w/^ [13] ■ Our 

solution is based on invariants of motion of the form [l^. 
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where Q = z /b -\- ga/uj^^ and Yi — bp — mbz + mb^gd/ujQ^ 
play the role of canonical variables. For Eq. ^ to be 
an invariant, the dimensionless functions b and a, re- 
spectively linked to the size az and center-of-mass po- 
sition Zcm of the cloud through az{t) = b{t)crz{0) and 
Zcmit) = —a{t)b{t)g/ujQz, must be solutions of 

d%/dt' + b{t)ujlit)^i,lMt)^ (3) 
d^a/dT^ + a{T) ^birf, (4) 
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FIG. 1: Trapping geometry (figure in the horizontal plane). 
Ultracold *^Rb atoms are trapped in an loffe-Pritchard-type 
magnetic trap created by current ig running through the 
three QUIC coils 1, 2, and 3. An additional pair of coils 
(a and b) produces an homogeneous field along y, which al- 
lows an independent tuning of the trap minimum field Bo via 
the current iBg- 



where T{t) — ujqz dt'/b^. The solutions of the time- 
dependent Schrodinger equation coincide with the sta- 
tionary states of the initial and final Hamiltonians H[t — 
0) and if (t/) if /(i = {0,t/}) oc H{t = {0,t/}) [11]. Thus 
we set d(0) = d(t/) = 6(0) = b{tf) = and a(0) = 1, 
a(tf ) — 7"^, 6(0) = 1, 6(t/) = 7. These latter terms im- 
ply that 6(0) — b{tf) — must hold as well, giving ten 
independent boundary conditions (BC). Our procedure 
to engineer uJz(t) is the following: (i) we use a polyno- 
mial ansatz for a(r) of the form a(T) = Sj=o '^ji'^/'^fV ^ 
for which ten coefficients are fixed by the I3C and the 
other can be arbitrarily chosen; (ii) we evaluate 6^(t) 
and thus 6[r(t)]; and (iii) using Eq. ^ we obtain the 
function ^^(i). Quite non- intuitively, the obtained solu- 
tion is valid for any magnitude of the linear term in the 
time-dependent Hamiltonian (as long as this linear term 
is not time-dependent itself). In the particular case of 
(7 = (no constant force term), however, a lower-order 
polynomial ansatz (fifth order) is sufficient [Tsj . 



III. EXPERIMENTAL PROCEDURE 

To experimentally investigate shortcuts to adiabatic- 
ity, we employ a sample of ultracold ^^Rb atoms held 
in a magnetic loffe-Pritchard trap. This popular type of 
trap is harmonic (for cold enough atoms) and anisotropic, 
with a typical ratio of 10 between the oscillation fre- 
quencies in the radial dimensions lJx^z and the axial one 
ujy (see Fig. [T} yielding the well-known cigar-shaped as- 
pect of the trapped cloud. For shallow traps, gravity 
significantly affects the potential in the vertical dimen- 



sion, yielding a displacement of the trap minimum —g/uj'^ 
compared to a tight trap. Our magnetic trap is of the 
quadrupole-Ioffe-configuration type (QUIC trap) intro- 
duced in Ref. [To'l, the three-coils setup sketched in Fig.[T] 
For sufficiently cold atoms (/csT -C uBo), the magnetic 
potential is harmonic of the form [2Q| 
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where fi/h « 1.4 MHz/G for our atoms in 
\F = 2,mF = +2). B' is the radial gradient of the mag- 
netic field while B" represents its curvature along y. Bq 
is the minimum of the magnetic field at the trap center, 
which can be adjusted using two independent parame- 
ters: the current ig running in the three QUIC coils, or 
the current iso in a pair of compensation coils providing 
a uniform field along y (see Fig.[T]). Since B" ^ B'^/Bq, 
the radial and axial angular frequencies are given by 
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These expressions show that we can, to some extent, ma- 
nipulate independently the radial and axial frequencies 
using iq and iso- 

Our initial sample is a small {N — 10^ atoms) and 
cold (To = 1.63 ^K) atomic cloud. The low temperature 
guarantees that the potential seen by the atoms remains 
harmonic even for large decompression factors. The small 
number of atoms is chosen to reduce the density and 
thus the elastic collision rate, responsible for the energy 
transfer between dimensions and thermalization. In the 
compressed trap with previously mentioned parameters, 
the typical time between two elastic collisions is « 28 ms, 
quite larger than the radial oscillation period of 4 ms. 

Prior to implementing a decompression sequence, we 
need to characterize the initial and final states. To this 
end, the position and size of the atomic cloud in three 
dimensions (3D) are measured using absorption imaging 
along two orthogonal directions. The trap frequencies 
are measured by slightly offsetting the trap center using 
compensation coils, then abruptly releasing it and mea- 
suring the cloud's center-of-mass motion as a function of 
time. For our fully compressed trap {ig — 27 A, iBg =0), 
we obtain vq^ — ujQx/2n = 228.1 Hz, voy = 22.2 Hz and 
i/Qz = 235.8 Hz. To measure the parameters of the fi- 
nal, decompressed state, we perform an "adiabatic-like" 
(i.e., slow tf = 6 s) decompression using linear ramps for 
the currents ig and iB„- In the following, we will refer 
to such ramps as "linear decompressions," although the 
resulting ujz{t) is not strictly linear [Eqs. ^ and ([7])]. 
The results presented in this article are obtained with 
a vertical decompression factor voz/vfz = 15, yielding 
final frequencies Vfx = 18.1 Hz, Vfy = 7.1 Hz, and 



Vfz 



15.7 Hz for the decompressed trap. In practice, 
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FIG. 2: Optimal trap frequency trajectory for a 35 ms verti- 
cal decompression. We plot (line) z^z(i) for a 35 ms vertical 
decompression from i/oz = 235.8 Hz to Vjz ~ 15.7 Hz, ob- 
tained with the invariant method (see text). The symbols 
correspond to measured values of the vertical trap frequency 
during the decompression process. 



this is achieved by decreasing iq from 27 to 3.6 A and 
increasing ibq from to 3 A. Since Vy is not affected by 
the increase of iso [see Eq. ([7])], the decompressed trap 
is much more isotropic {yf{x,z}/^fy ~ 2) than the com- 
pressed one {i^o{x,z}/i^Oy ~ 12). 

We illustrate the efficiency of our shortcut method by 
realizing a fast {tf = 35 ms) trap decompression opti- 
mized for the vertical dimension 2, where gravity strongly 
affects the cloud's motion. The employed solution z^z(i) is 
shown in Fig. [5] (line, note the vertical log scale). Because 
of the finite time response of the trap electronic circuit, 
the measured trap field profile is different from the com- 
puted one. We thus monitored Vz by interrupting the 
sequence at different times, and adjusted the compensa- 
tion field to obtain a measured i^zit) (symbols in Fig. [5]) 
close to the theoretical one (deviation < 5%). The un- 
certainty on the experimental values is ±2%. 



IV. RESULTS 

Figure [3] shows the result of the shortcut decompres- 
sion using the trajectory of Fig. [2] We plot in Fig. ^a.) 
the time evolution of the cloud's center-of-mass position 
Zcm once the decompression sequence is completed, and 
in Fig. [njb) that of the cloud's size (Jz- These data cor- 
respond to averages over three successive images, taken 
after a 6 ms time of fiight. The open circles correspond to 
an abrupt jump from vqz to Ufz (in practice, the effective 
decompression time is « O.I ms). The solid circles are ob- 
tained with a (non-optimal) linear decompression of du- 
ration tf = 35 ms, the stars with the shortcut trajectory, 
and the squares with a quasi- adiabatic linear decompres- 
sion in 6 s. In every instance, we observe in Fig. [SJa) 
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FIG. 3: (Color online) Vertical trap decompression: compar- 
ison between different schemes. We report in (a) and (b), 
respectively, the cloud's vertical center-of-mass position Zcm 
and size versus time after decompression, for four different 
sequences. Open circles (green): abrupt decompression; solid 
circles (black): linear decompression in 35 ms; stars (red): 
shortcut decompression in 35 ms; squares (blue): linear de- 
compression in 6 s. 



the expected sinusoidal oscillations of Zcm at the decom- 
pressed trap frequency Vfz ~ 15.7 Hz (dipole mode) and 
of amplitude Azcm- As can be seen, the shortcut decom- 
pression yields a strong reduction of Azcm when com- 
pared to the abrupt and 35 ms linear decompressions, 
by a factor 9 and 7.2, respectively. However, the resid- 
ual center-of-mass oscillations after the shortcut sequence 
are still sizable, a factor of 5 larger than that observed 
for the 6-s-long linear decompression. We attribute these 
residual oscillations to imperfections of the experiments 
which are discussed at the end of the paper. From the 
amplitude Azcm we can infer the excess energy commu- 
nicated to the cloud in the form of the dipole excitation 
Edip = l/2mu}'j^Azcm'^ . We also observe in Fig. |31^b) 
oscillations of the cloud's size CTz at twice the frequency 
of the decompressed trap (breathing mode) . The theory 
predicts such non-sinusoidal periodic oscillations, whose 
expression can be derived analytically ^21j . In the exper- 
iment, the measurement of cr^ is less accurate than that 
of Zcm , because of the limited spatial resolution and noise 
{uz ~ 40/im), and we cannot fit the measured oscillations 
to the model in every instance. We thus quantify the am- 
plitude of the breathing mode by using the standard de- 
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viation Actz of C7z{t) after decompression. We observe a 
reduction of Acr^ when we use the shortcut trajectory, by 
a factor 7 and 3 when compared to the abrupt and hnear 
decompressions respectively. The residual Aa^ is again 
a factor of 5 above that of the 6-s-long linear decompres- 
sion. The excess energy stored in the breathing mode is 
Ebreath ~ 2mu>'j^Aaz^ ■ Note that we overestimate Acr^ 
(and thus Ebreath) because of our 6 ms time of flight. 
The total excess energy imparted to the system during 
the decompression is then -Eexc = -E-dip + E{,j,eath ■ Quite 
obviously from the vertical scales in Fig. |31 we always 
have Edip ^ Ebreath- Using the previous expression, we 
find excess energies of 54, 35, 0.7, and 0.02 /zK for the 
abrupt, 35 ms linear, 35 ms shortcut, and 6 s linear de- 
compressions, respectively. For the latter, we measured a 
final temperature Tf — 0.13 ^K. Since the initial temper- 
ature is Tq = 1.63 fiK, the cooling factor is 12.5, quite 
close to the expected v^zlvfz = 15 value for a purely 
adiabatic transition. 

To provide the reader with a better feeling of the time 
scales involved in the trap decompression, we compare in 
Fig. m our shortcut results with those of linear decom- 
pressions with various durations (full circles for Azcm, 
open circles for Aa^. All the amplitudes in this fig- 
ure are normalized to those corresponding to an abrupt 
decompression [tf = 0.1 ms). The stars correspond to 
three shortcut experiments. Two experiments were per- 
formed along the vertical: the 35 ms one depicted on 
Figs. [2] and [31 and a 100-ms-long one. Another, 20-ms- 
long shortcut decompression was also performed along x 
(no grav ity) , using the fifth-order polynomial ansatz of 
Ref. [13|. The solid stars correspond to oscillation am- 
plitudes of Zcm while the open stars stand for Aa^ ■ The 
linear decompression data allow us to estimate a quanti- 
tative criterion for adiabaticity instead of the usual qual- 
itative criterion t/ 3> For instance, we can set as 
a criterion that the excess energy should be of the or- 
der or smaller than the thermal energy associated with 
Tf = To/72: ~ l/2mu;j^Az,^^ < ksTf- This 
condition yields tf > 3.3 s for linear ramps. Since our 
shortcut sequence in 100 ms also satisfies this condition 
it can be considered adiabatic using this criterion, with a 
reduction of the necessary decompression time by a fac- 
tor of 33. The 35 ms shortcut decompression does not 
meet the criterion above but still realizes a gain on the 
transition time of a factor of 37 when compared to a lin- 
ear ramp. The 20-ms-long decompression along x excites 
only the breathing mode (no trap displacement), whose 
residual amplitude is one order of magnitude lower than 
for the abrupt transition, and a factor of 2 above that of 
the 6 s linear ramp. 

As stressed in the theoretical part of this article, only 
the final state is identical to that obtained through an 
adiabatic sequence. Indeed, we performed an experiment 
where we interrupted the frequency trajectory of Fig. [2] 
after 10 ms (as pointed out by the arrow in the figure). 
Despite the fact that 94% of the frequency difference 
vqz — Vfz has been covered at t = 10 ms, we observe large 
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FIG. 4: Summary of faster-than-adiabatic decompression re- 
sults. We plot the amplitudes of center-of-mass (filled circles) 
and cloud's size (open circles) oscillations along the vertical 
direction after linear decompressions of various durations tf. 
All amplitudes are scaled to that of the abrupt decompression 
{tf = 0.1 ms). The stars correspond to our shortcut decom- 
pression experiments in 20, 35 and 100 ms (filled symbols: 
center of mass, open symbols: size). 



center-of-mass and cloud size oscillations, respectively a 
factor of 7 and 4 larger than those observed when the 
entire 35 ms sequence is completed. Thus, the last 25 ms 
of the frequency trajectory in Fig. [5] are of paramount 
importance for reaching the optimal final state. 

We now discuss experimental imperfections which 
might be responsible for the residual oscillations observed 
in Fig. [3l The first possible cause is a mismatch between 
the theoretical frequency trajectory and the experimental 
one. As shown in Fig.[2l we did our best to maintain this 
mismatch below 5% for selected time values of the tra- 
jectory, but we cannot guarantee that this holds for the 
whole sequence. In particular, as discussed previously, 
the last part of the trajectory where the frequencies are 
small and thus the relative measurement error large is po- 
tentially more critical. Probably most importantly, our 
trap can be considered harmonic only for small atomic 
displacements from the trap center. During the short- 
cut decompression, the trap center shifts vertically by 
1 mm and the atoms follow a complex dynamics that 
brings them quite far from the trap center (« 300/im). 
Deviations from harmonicity may thus play an important 
role in our experiment 11], limiting the performances of 
our shortcut decompression. 

V. CONCLUSION 

In conclusion, we presented in this article the first ex- 
perimental realization of the faster-than-adiabatic dis- 
placement and cooling of an ensemble of magnetically 
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trapped ultracold atoms using an optimal decompression 
sequence based on invariants of motion. Using this for- 
malism, we derived optimal trap frequency trajectories 
in the case of a time-dependent harmonic potential plus 
a time-independent linear term accounting for gravity. 
Our solution also applies to the simpler case of a purely 
harmonic potential such as that treated in Ref. '13]. We 
demonstrated the validity of our scheme by applying a 
fast (35 ms) 15-fold frequency decompression to the trap 
in the vertical dimension, yielding a residual center-of- 
mass oscillation of the cloud equivalent to that of 1.3-s- 
long linear decompression (a reduction by a factor of 37). 
As a future prospect, one could apply this technique to 
more isotropic traps (such as crossed dipole traps) to 
obtain a faster and efficient cooling in 3D and produce 
very low temperatures. Optimal trajectories could also 
be searched for in other situations such as the moving 
quadrupole magnetic traps often used to transport cold 



atoms [7[. This method can also be readily applied to 
a Tonks gas [2l|, and to Bose-Einstein condensates with 
some restrictions on the dimensionality due to the scal- 
ing of the interaction term [1^ . More generally, these 
optimal faster-than-adiabatic schemes could be adapted 
to many areas of physics where time-dependent Hamil- 
tonians are employed. 
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